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Abstract 

The asymptotic behaviour of a closed BCMP network, with n 
queues and m n clients, is analyzed when n and m n become simul- 
taneously large. Our method relies on Berry-Esseen type approxima- 
tions coming in the Central Limit Theorem. We construct critical se- 
quences m° , which are necessary and sufficient to distinguish between 
saturated and non-saturated regimes for the network. Several applica- 
tions of these results are presented. It is shown that some queues can 
act as bottlenecks, limiting thus the global efficiency of the system. 

1 Introduction 

In many applications (telecommunications, transportation, etc.), it is desir- 
able to understand the behaviour and performance of stochastic networks as 
their size increases. From an engineering point of view, the problem can be 
roughly formulated as follows: 

Consider a closed network with n nodes and exactly m n customers 
circulating inside. Find a function f , such that m = f(n) yields 
an interesting performance of the system as n increases. 

In this study, we start from the so-called product-form networks, which 
play an important role in quantitative analysis of systems. Although the 
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equilibrium state probabilities have then a simple expression (see for example 
Kelly [1]), non-trivial problems remain, due to an intrinsic combinatorial ex- 
plosion in the formulas, especially in those involving the famous normalizing 
constant. To circumvent these drawbacks, the idea is to compute asymptotic 
expansions of the characteristic values of the network, when m and n both 
tend to infinity. 

This approach has been considered by Knessl and Tier [5], Kogan and 
Birman [EHHd] and Malyshev and Yakovlev |10J . However, it relies on purely 
analytical tools, which are difficult to use in a more general setting and, in 
our opinion, do not really give a structural explanation of the phenomena 
involved. 

The method proposed hereafter has direct connections with the Central 
Limit Theorem: instead of representing the values of interest as complex inte- 
grals, we express them in terms of distributions of scaled sums of independent 
random variables. Besides giving a clear interpretation of the computations, 
this allows to handle directly the general case of single-chain closed networks. 
We show by construction the existence of critical sequences m a n in the follow- 
ing sense: the network saturates if, and only if, m n 3> m° . These results can 
also be interpreted as insensitivity properties: as the number of stations n 
and the number of customers m n go to infinity, the network is shown to be 
equivalent to an open network of n independent queues (having a total mean 
number of customers m n ), in the sense that both systems have asymptotically 
the same finite-dimensional distributions. 

The paper is organized as follows. The model is introduced in Section [21 
together with a presentation of the method. In Section [3j asymptotics of the 
marginal distribution of the queue lengths are given under normal conditions 
and also when some queues become overloaded. Section H] unifies the results 
and contains the main theorems about scaling. Section \5\ and El are devoted to 
concrete applications of these results, in particular to service vehicle networks 
(like the Praxitele project, now developed at INRIA). Section [7| contains some 
conclusive remarks. Most of the technical proofs are postponed in Appendix. 

2 Mathematical model and view of the main 
results 

Consider a closed BCMP network C n with n queues and m n clients. The 
number of clients at queue k at steady state is a random variable Qk, n - The 
service rate at queue k when there are qu customers is fik,n{Qk)- The routing 
probability from queue k to queue t is Pk,e, n and P n denotes the transition 
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matrix supposed to be ergodic, with invariant measure ir n = (7r l n , . . . , 7T n)fl ), 
defined by: 

n n P n = 7T„ and 7i hn H h 7r nin = 1. (2.1) 

Then it is known that, for any qi, . . . , q n > such that q± + ■ • • + q n = m n , 

P.«U = *,..., - *) - Z-l„ n ^ (1) . t ; tnte) . (2.2) 
with the normalizing condition 



Z 



7C Qk 



m,n 



TT (2.3) 



(JiH hgn=m h 



It is worth noting that our analysis applies to any network which has a 
product form equilibrium distribution like (12. 2p . It includes for example, as 
soon as the matrix P n is reversible, all systems having transition rates of 
the form Pk,e,n&k,n(qk) (3e,n(qe) , in which case finite capacity situations can be 
covered (e.g. (3i tn (qi) = 1{^<^})- See Serfozo [TT] for further examples. 

To avoid hiding global results with tedious technicalities, we suppose 
throughout the study that, for all n, C n contains at least one queue which, 
taken in isolation, can be saturated with a finite input flow (e.g. a MjMj cj oo 
queue) . 

The overall presentation requires definitions and an intermediate lemma, 
given in Section 12.11 The informal presentation of the central results appears 
in Section |2~2"1 



2.1 Preliminaries 

Define, for each k, 1 < k < n, the generating function 



fk,n( Z ) = Yl 



z q 



Q //fc,n(l) • • • VkM 

Note that for each n, fk, n has a singularity at finite distance for at least 
one 1 < k < n. 

To the original closed network C n , we let correspond a new system O n {\) 
which is open and consists of n parallel queues, with service rates (ik,n(x) and 
arrival intensity \~Kk,n at queue k, where the choice of A will be made more 
precise later. The queue length Xfc jn (A) of the fc-th queue of O n (X) has a 
distribution given by 

P(X k>n (X) = x) - 



,n) f^k,n 

(1)" 
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and X l n (A), . . . , X n n (A) are independent variables. We assume that X k n {\) 
has some finite moments of order r > 2 and introduce the following notation: 

m kjn (X) = EX k , n (X), S n (X) = Ylk=iXk,n(X), 

C(A) = E|X M (A)-m*,„(A)r, )9«(A) = ELi^(A), 

= E[X M (A)-m M ,(A)] 3 , 0«>(A) = YT^J&W- 

Let <Pk,n(9 ; A) be the characteristic function of X fc „(A) — m fc „(A). Then, 
for any real 9, 

<p k>n {9 ; A) = Ee K**,»(A)-"*,»W)« = ^^^ -^W, (2.4) 

and 

; A) d = f Ee H^(A)-E5 n (A) )e = ^ . A) . . . . A) (2 5) 

The reason why O n (A) has been introduced is that the main performance 
characteristics of the network C n can be expressed simply in terms of the 
distribution of Xl jH (A), . . . , X Hjn (A): 

Lemma 2.1 (i) For any choice ofm n , there exists a unique X n such that 
ES n (X n ) = E[X 1>n (X n ) + ■■■ + X njn (X n )) = m n . (2.6) 

From now on, unless otherwise stated, all quantities will pertain to the 
network O n (X n ) and X n will be omitted. 

(ii) Equations 112.2]) and 112.8]) can be rewritten as 

1 n 

P(<9l,n = <?!,■■•, Qn,n = q n ) = (q F(^fc,„ = q h ). (2.7) 

(m) For any i > and qi, . . . ,qe > ; the joint distribution of the number 
of customers in the queues 1, . . . ,£ of C n is 

F(Qi,n = qi,---,Qe,n = qe) (2.8) 
H S n - J2k=i X k,n = m n - YjI=i Ik) 



F(S n = m n 
j • • • ) ^e,n — 
and, consequently, EQe, n = E[Xg jn \S n = m n ]. 



l[HXk,n 



Qk) 



k=l 

F(X hn = q x , . . . , Xt, n = qe\S n = m n ), 
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[iv] For any 1 < I < n, 

ivn gO 5 ^ -X e + X e = m n ) 

Efe = ^ P(£ n = m n ) ' (2 ' 9) 

where X^ n is an integer-valued r.v., independent from everything else 
and having distribution 

F{X i>n = x) = . 

m£ tn 

Note that A n can be obtained as the unique solution of the equation 



m n = 

k 



v fc '7 ra ; fc -" v ; fc '" n; . (2.10) 



While this equation is in general impossible to solve explicitly, A n can be 
computed numerically using classical methods. 

Proof A straightforward computation yields, for all 1 < k < n, 

dm k ,nW a k,nW 



d\ A 



> o. (2.11; 



The mean number of clients in O n (X) is thus a strictly increasing function 
of A, which equals zero when A = and goes to infinity with A. This proves 
the first assertion of the lemma. 
Define 



Y, 



7 \ m 

def ^nin^n^n 



n i tti 



Wk=l fk,n\^k,n^r 

Then reads 

1 - 

Pn(Ql,n = Ql, ■ ■ -,Qn,n = Qn) = JJ H X k,n = Ik), 

which yields (12 .7p . since 



n,n " "n 



qi H h<?n=«ln k = l 
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Equation (12. 9p and the first part of (12. 8p are derived similarly. For the 
second part of (I2.8p . we simply note that 

i 

H S n ~ J2k=l X k,n = m n - Yj[=1 9k) Yl H X k,n = Qk) 

k=l 

I 

= F(S n = m n \X hn = qi, . . .,X t>n = q t ) JJp(X fci „ = q k ) 

k=l 

= F(X hn = qi, . . . , X £y n = qi\S n = m n ) F(S n = m n ). 



2.2 Informal description of the method 

Most of the derivations obtained in the paper are based on the various repre- 
sentations given in Lemma 12. 11 Whereas the studies O El E EE] use mainly 
saddle-point methods, our approach relies on direct limit theorems for the 
distribution of S n . 

For example, assume that S n — m n satisfies a local limit theorem such as: 

Under "suitable" conditions, there exists a distribution with den- 
sity f and a sequence a n such that, for any integer x, 

lim a n F(S n -m n = x) - f(-) = 0. (2.12) 

n-s-oo ^O n / 

Then Lemma 12.11 will yield 

n 

HQhn = Ql, ■ ■ ■ , Qn,n = Qn) ~ J^j JJ H X k,n = qk) , 

and, for any finite £, 

i 

HQi,n = qi,>>>,Qi,n = qi) ~ Y[H x k, n = Qk)- 

k=l 

This amounts to say that the joint distribution of any finite number of 
queues in the BCMP network C n is, at steady state, asymptotically equivalent 
to the product distribution of the corresponding queues in the system O n . 

It is at this moment important to emphasize that we do not require any 
"smooth" limiting behaviour for O n , which is somehow an instrumental net- 
work, computationally easier to evaluate. 
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To prove local limit theorems like ( I2.12p . it is necessary to investigate 
carefully the behaviour of the variables X^^. In particular, since ~ES n = 
m n < oo, all queues in O n are ergodic, which reads, for any 1 < k < n, 



An7Tfc, n < Pk,n < OO 



or, equivalently, 
where typically 



p° n = \ n max ^ < 1, (2.13) 



Three main situations have been analyzed: 

i) p° n is bounded away from 1: then S n /a n satisfies a local Central Limit 
Theorem and tends in distribution to a normal law (see Theorem I 



(ii) p° n — y 1 and the supremum in i2.13\) is attained for a finite number of 
queues: then the network subdivides into two subsets, the "saturated" 
queues and the rest of the network. As shown in Theorem I4.3[ under 
mild regularity assumptions, there exists a sequence a n such that S n /a n 
tends to a gamma law; 

(Hi) p° — y 1 and the supremum in 112. lty) is attained for an unbounded num- 
ber of queues: S n /a n again tends in distribution to a normal law (see 
Theorem 14.4)1 . 

In fact, Theorems 14.2) I4T31 and H~4l quoted above are general, in the sense 
that they provide a construction of efficient scalings in terms of m n , the 
number of customers: the existence of critical sequences m° for the network 
C n is shown by explicit construction. Under reasonable assumptions, these 
sequences are necessary and sufficient to discriminate between saturated and 
non saturated regimes. This is similar to phase transition phenomena ob- 
served in [TO], where it was assumed that m n /n — > A > (see Section IfTTj) . 
Clearly, for a non-saturated regime to exist as n — > oo, it is necessary to have 
fn n = 0(n); this condition is not sufficient (see Section 1672]) . 

Condition (12. 13)) can be used to determine an upper bound for A n and to 
exhibit queues which act as bottlenecks in the network C n (see Section H]). 

Remark Rather than simple limit theorems, the results in Sections [3] and S] 
are given in terms of asymptotic expansions, using the operators O and ft 
defined as follows: 

a(n) = 0(6(77)), iff 3iT>0, V77, \a(r])\ < K\b(r])\, 
0(77) = 0(6(77)), ^ a (v) = 0(b(r})) and b(r]) = 0(0(77)), 
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where 77 is some unspecified argument. Unless otherwise stated, all these 
bounds are uniform with respect to n and all queue indexes. 

3 Local limit theorems and asymptotic ex- 
pansions 

In this section, we compute estimates of several performance measures of 
C n by means of local limit theorems on sums of independent random vari- 
ables. The two series of results presented here are of somewhat different na- 
ture: whereas the conditions of Proposition 13.11 depend on moments, Propo- 
sition 13.31 relies on analytic properties of the generating function of some 
queues. 

3.1 Normal traffic case 

When the queues are not saturated (in a sense made more precise in Proposi- 
tion [3J]), it is possible to prove local Central Limit Theorems, relying more 
exactly on Berry-Esseen type expansions (see for instance Feller [3]). 
Define 7I n from X^^ n as in Lemma IA.1I of the appendix, and let 



Proposition 3.1 (i) Let, for any < r < 1 such that P n 2+r) exists, 



Let 7„<5 n — > 00 as n — >■ 00. Then, for any integer x, the following 
approximation holds uniformly in x: 



2 dcf 2 I ,2^2 
In = 7l,n + • ■ ■ + 7n,n < °n 



(V 
n 



del 



<7„P(S, 




1 



e 




(3.1) 



n 




°(!p) +0 ( 




n 



exp 




) 



) 



(m) Let, for any < r < 1 such that /3^ 3+r) exists, 



5, 



n 



dcf 
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Let 7„5 n — > oo as n — > oo. Then, for any integer x, the following 
approximation holds uniformly in x: 



o-n F(S n - m r 



1 



x) — 



2vr 
O 



i + gg> /x 3 3 X 



(3+r) 



a 



3+r 



?i " n 



(3.2) 



ll 8 n 



exp 



Proof See Appendix IA.2I ■ 

The main assumption of the previous proposition is classical, since it is 
nothing else but Lyapounov's condition, popular in the Central Limit Prob- 
lem: 

fl(2 + r) 

for some r > 0, lim = 0. (3.3) 



n— >oo (J 



This condition yields in particular (see e.g. Loeve [9]) 



lim max — 

n— >oo l<fc<n (J„ 



0. 



(3.4) 



which in turn implies the uniform asymptotic negligibility of the X^s. Note 
that it would be possible by truncation methods to prove similar results 
without requiring the existence of moments. 

We are now in a position to present some basic estimates when the size 
of the network increases. 



Theorem 3.2 Let r be a real number such that < r < 2. Assume that 
o n = 0(-y n ), that (3 n 2+r) exists and f3 n 2+r) / a^ +r — >■ as n — )■ oo. Then the 
following asymptotic expansions hold. 



P(Ql,n = 91, • • - ,<5n,n = <?n) 



27T<7 n JJ P(X fc , n = q k ) 



1 + 



(2 + r) 



fc=l 



a 



2+r 



(3.5) 



[ii) For any finite £, if ELi - tyl/ffn -» 0, 



P(Ql,n = 91, • • • , <?£,n = ft) = n P ( Xfc -« = ^ 



fc=l 



1 + 0(e lt . 



(3.6) 
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£ 1,< 



2+r 



+ l{r>l} 



(E,-=i m j,n-qj)fi. 



(3) 

n 



^ra) For any j ; 



(3.7) 



^2,n — — TTTT H ~ + 



(2+r) 



/? (3) / a- 



0" 



2+r 



m i,n<T l n +r 



Proof Equation (13. 5p is a simple application of Proposition 13. II to (12. 7p . 
To prove (13. 6p from (12.81) when r < 1, we simply write 

F( S n ~ J2k=l X k,n = m n - Yll=l Qk) 



F(S n = m n ) 



-L 9 



/P(2+r) (sZj = l m 3,n-tj) 

l + 0{^— + e ^ 



2+r 



- 1 



and use the relation |e " 2 / 2 — 1| < u 2 . When 1 < r < 2, it suffices to take 
into account the inequality 

(n 3 - 3 M )e-" 2/2 = 0(u). 

Relation (13.71) is also derived from (12.81) . ■ 



3.2 Heavy traffic case 

We proceed now to analyze the behavior of the network C n when some queues 
saturate, as n — > oo. This, in particular, implies that the Lyapounov condi- 
tion ( 13 .3p is no more valid. In fact, after a suitable normalization, S n —m n will 
be shown to converge in distribution to a random variable having a gamma 
distribution, under the broad assumption that the first singularities of the 
relevant generating functions are algebraic. 

Let, for some p° n G [0, 1[ (to be specified in Section HJ), 

l_ p o el 6P 

and assume 
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Assumption Al There exists a set J 7 ® of "saturable" queues, such that, for 
all k G there exist a real number ^f.,n ana " a function ipk,n{6) satisfying 
the relation 

Moreover, ip' kn {0) = 0(1), uniformly in k and n, and there exists a con- 
stant £ max such that 

1 — <,k,n *^ £max ^ OO- 

Clearly, the term cut'" (6) coming in the definition of (#) emphasizes 
the fact that the generating function fk, n (z) pertaining to queue fceTJ has 
its first singularity which is algebraic of order £ fc ,n- If, in addition, p° n — > 1 
as n — > oo, the working conditions of the system ensure all queues in J 7 " 
saturate so that, in particular, E,Xk t n ~ £fc,n«n, where 

n° 

dof P n 
Ctn - 



While this assumption covers a wide range of known queues, it is clear 
that other types of singularities could be handled via the same method. 
Let 

and define the total characteristic function of the queues in T n \ J 7 " by 

= n 

Let r be a real number, < r < 1. Hereafter, a n , j3£ +r) , % and 5 n will 
denote quantities having the same meaning as in Proposition 13. 1[ but related 

tO (fi n {9). 

The counterpart of Theorem l3.2l now reads, in the case of heavy operating 
conditions: 

Proposition 3.3 Let p° n — >■ 1. 7/£ n is bounded, a n /a n — > and 5 n 7„ — > oo 
as n —t- oo, then the following estimate holds: 

a n F(S n -m n = x)- - (3.8) 

Of (*»Y + -1 + ^l(^) 2+r + Pn +T) f&n^- 1 



+0 ( 



e s 
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Proof See Appendix IA.2I 



The estimates of Proposition 13.31 allow to establish the main result of this 
section. 

Theorem 3.4 Let a n /a n -» 0, j3£ +r) /cfl +r -»■ and a n = 0(%) as n -> oo. 
Then the following expansions hold when £ n is uniformly bounded: 



[i) for any q h ...,q n >0, 



HQl,n = Ql, ■ ■ -,Qn,n = <?n) 

««r(£ 



J JJ H X k,n = q k ) 



k=l 



l + 0(e r 



with 



a,, 



+ — + 



+ 



cr 2+r \a 

u n 



(3.9) 



[ii) for any finite I, such that fl [1, £] = 0, 

HQi,n = qu---,Qe,n = qe) (3.10) 

= fl = 9k) [l + 0(e n ) + O (^ =1 mKn " 9fc ' 

fc=i 



a, 



^m) /or any j & 



EQj,n = EXj 



3,n 



1 + 0(e„) 



(3.11) 



» /or any j <G J£ ; 



EQj,n = EX jjn 1 + 0(£ n 



(3.12) 



Proof The proof of (i) -(Hi) follows essentially along the same lines as for 
Theorem I3.2[ while (iv) depends on Equation ( 12. 9 j) of Lemma [2.11 ■ 
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4 Scaling 



As said in the introduction, this section provides guidelines for using the 
above technical results in two ways. 

• Quantitative estimates for the error terms (w.r.t. some limiting distribu- 
tion), explicitly obtained from the original data (e.g. the total number 
of customers m n ). 

• Qualitative understanding of the "critical" values for m n which, in some 
sense, induce phase transitions of interest. 

The queues are partitioned as follows: 



Tn = \0<k<n : lim v7/i fcin (l) • • • H k ,n{q) < oo L 



In = \0<k<n : lim ^/// fc)n (l) • • • fi k>n (q) = oo \. 

From the general discussion at the beginning of Section |2j T n is never 
empty. Let also 



dcf 



lim f//i fc)Tl (l) • • -fik,n(Q), tfkeFn 

Pk,n — \ q— >oo 

Mfc,n(l), like!, 



n ■ 

dcf ^n^k,n ■. n dcf . Pk,n n dcf ^n, 

pk,n = , K = mm , p n = maxp fci „ = — . 

Pk,n keJ- n 7Tk t n k<=J- n A n 

We shall also need the following subset of T n : 

K = {keFn : p k)U = p° n }. 

Note that the definitions of p k ,n and p° n are consistent with the discussion 
which lead to (12.13p . Moreover, in most practical cases, Pk,n(q) — > Pk,n as 
q — > oo, provided that this limit exists and is finite. 

To avoid uninteresting technicalities, it will be convenient to introduce 
Assumptions IA2I and IA3[ but it should be pointed out that the results of 
Section [3] are valid in a more general setting. Simple conditions ensuring IA1I 
and [A3 are discussed in Section |5j 

Assumption A2 The following limit holds: 

lim max 



n^-oo l<k<n — 1 r " " " ~r — 
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Assumption IA2I is somehow unavoidable to obtain a meaningful asymp- 
totic behaviour of the network. It says that it is possible to let m n — > oo 
as n — > oo, without saturating the network and, under the forthcoming As- 
sumption [AS it amounts to Lyapounov's condition (13.31) . Note that, when 
fik,n = ^(1) uniformly in k and n, IA2I is simply equivalent to 

lim max 7Tk, n = 0. 

n— >oo l<k<n ' 

Assumption A3 (i) For any real A < 1 and any integer r < A, and for 
any k G T n such that pk, n < A, 

m k ,n = M(Pk,n), = ^(Pk,n), Tfc.n = ^(pk,n) (4.1) 

uniformly in k and n. 
(ii) d^.l[) also holds for all k e I n . 



The derivation of the most general results of the section is done in Lemma 
14.11 and Theorem 14.21 Further insight, under some additional assumptions, 
is presented in Theorems 14.31 and 14.41 

Definition A sequence m° n is said to be weakly critical for C n if, for any 
< t < 1, 



n 



exists and lim g(t) be either 1 or oo. 
If, in addition, the relation 



lim lim mn{tK) = lim iinT 



holds, then the sequence is said to be strongly critical for C n . 

Before seeing how such critical sequences can be used, the next lemma 
proves their existence. 

Lemma 4.1 Under assumption \AZ\ a convenient weakly critical sequence 
for C n is, for some fixed < u < 1, 

m° n (u) = h u m n (u\° n ), (4.3) 

where h u is correctly chosen. 
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Proof Choose (t,u) e]0, l[x]0, 1[. From[A3 



l^k,n 



k=l 

and the application t !->■ m n (tA°)/m n (?/A°) is increasing and locally bounded. 
Therefore, 

* ,,x def ^n(^A°) 

# u (t) = hm — — 

n->oo m n (wA°J 

exists and is increasing. To conclude the proof, take 



h, 



lim g u {t), if the limit is finite, 
1, otherwise. 



It is interesting to note that, if the above limit is finite for some u, it is 
finite for all u e]0, 1[. The proof of the lemma is concluded. ■ 



In fact, as shown in Theorem I4.2[ any critical sequence m° acts as a 
threshold parameter for m n . Under IA2I and IA3| which are satisfied by a 
wide variety of networks, we provide a nearly complete classification in terms 
of necessary and sufficient scaling. It is worth to emphasize that any m° n 
chosen from ( 14. 3 p has a pseudo- explicit form, given in terms of the data of 
the original network. 

The second step is to enumerate in a consistent way the desirable proper- 
ties of the distribution of Qi ;Tl , . . . , Q n ,n- for some finite j and some unspeci- 
fied e n , such that e n — > as n — > oo, we have 



EQj, r . 



EX. 



[l + 0(e 



(4.4) 



HQi,n = qi,---,Qj,n = qj) = Y[H X k,n = qk)h + 0{e n )\, (4.5) 



fc=i 



and also, when Theorem 13 . 2 1 [resp . Theorem 13 .4] holds, the following equation 
2D [resp. flU}]: 



HQl.n = Ql, ■ ■ -,Qn,n = ?n) 

n 

= V^rtJn JJp(X fc , n = q k )\l + 0(e n ) 



k=l 



P(Ql,n = 91, • • -,Qn,n = Qn 



^Unx k ,n = qk) 



k=l 



l + 0(e r 



(4.6) 



(4.7) 
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Theorem 4.2 Let \A2\ and \A3\ hold and m° n be a weakly critical sequence for 
C n , with the associated function git). 

Assume first that lim^i- g(t) = 1. Then the following classification 
holds: 



hm < 1, 



then \4-4\) > H4-5j ) an d \4-ty hold with e n = l/m n . In particular ~EQk, 
is bounded, uniformly in k and n. 



hm > 1, 



then, for any sequence of queues k n in J^°, we have lim KQk n ,n — oo. 
(Hi) If m° n is a strongly critical sequence and 

lim — - > 1, 

then, for any sequence of queues k n in T^, we have lim W,Qk„,n = oo. 

In the situation lim^i_ g(t) = oo, the same results hold, just replacing 
"< 1" (resp. "> V ) in the r.h.s. of the inequalities by "< oo" (resp. "= oo"). 

Proof To prove (i) , note that m n = m n (X n ) = m n (p° A°). Since m n (tA°) 
is increasing in t, this implies that, when lim^oop" = 1, we have also 
lim^oo m n /m° n > 1. Therefore, in case (i) , there exists r < 1 such that 
fP n < r for any n£l Using IA31 we can estimate all error terms coming in 
Theorem 13.21 and the result is proved. 

Similarly in case (it) [resp. (Hi) ], we have necessarily lim^oo — 1 [resp. 
lim^oo p° n — 1] , and the result follows from the monotonicity of the function 
t m- m kn , n (t\° n ). 

The case lim t _>.i_ g(t) = oo is handled with the same method. ■ 



Direct applications of Theorem 14. 21 are proposed farther on in sections 16.1 1 
and 16.21 

In order to get finer results, the next assumption ensures that the queues 
not belonging to J-° stay uniformly away from saturation conditions. 



Assumption A4 There exists a constant A < 1 such that, 

''Hk 



A° ^ < A, for all k e J*„ \ J*, (4.8) 
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In order to properly reformulate the results of Section [31 let us define 

^n(A) = ^m fc , n (A), (4.9) 

K = ™ n (K)- (4.io) 



Using (12. lip , it is not difficult to see that m° n defined (I4.10p is a strongly 
critical sequence for C n under lAlj. \A2\ IA3I and IA4I Therefore, all results of 
Theorem 14.21 hold, as well as the following: 

Theorem 4.3 Let\Al\ \A2\ IA3I and IA4I hold. If £ n is uniformly bounded, 
then the following results hold: 

(i) If there exists 9 n > 0, such that, for all n <E N, 



rn 



and lim^oo 0%m n = oo, then \4-4\) > an d M-°D hold with 



dc 



f 1 1 



m„ ml 9i 



except when a queue in J 7 " is concerned, in which case \4-4\ ) an d 
hold with 

1 

d l m n ' 

[ii) If there exists 9 n > 0, such that, for all n6N, 

-i>i + e n , 

and lim^^oo 6* n m° = lim n ^oo O^rh^ = oo, then \4-4\) an d M-ty hold with 

1 1 1 .- 1 Sn-l 

def J- 

~ n 



m°J 2 



n n w n" u n \J llb n n n 

Moreover, if in Equation ( |^.<5| ), [IjiJnJ 7 " = 0, then the latter also holds, 
with e n having the above value. 

Proof To prove (i) , note that when m n < (1 — 9 n )m° n , 

m n (\ n ) < m n (X n ) < (1 - n )m n (X° n ). 
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Moreover, using \A3\ IA4I and ( 12. lip . Taylor's formula yields, for some 
A G]A n , A° [, 

m° n -m n (\ n ) = m n (A°) - rh n (\ n ) 



(A° - A n ) 



A 



which implies 
Hence, 



1-^^ = 0(1-/*) >*„. 



or 1 



1 - P°„ 

and, using m n (A n ) = f2(/3^ r) ) = Q(a%), a direct but tedious computation 
shows that Theorem 13.21 applies with appropriate error terms. 
Let us now prove assertion (ii) . It follows from 

Pi 

that p° — > 1 and 



m n - rh n (\ n ) = ( ^"^"o ) > ^n^n °°, 



^ = fi(m„(A n )(l-p^ 



(m n - m n (X n )) 2 J \6l[rh° n } 2 J V^m 

Thus, Theorem 13.41 applies and (ii) is proved. ■ 

It remains to state what happens when £ n — > oo as n — > oo. As shown 
below, this behaviour does not depend on the saturation of the queues in J 7 ". 

Theorem 4.4 Let £ n — > oo as n — >■ oo. Let a/so IA1|. [A2L IA3I and I A4I /ioW. 

Then, under the uniformity assumption 

C = ((]%f)i) ' for dl k G ^ ^ 

the results h4-4\ ), M-5\ ) and are again i>a/id, with 



def 



1 - Pn) m n 
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Proof The statement relies on Theorem 13.2^ taking r = 2. First, from 
classical weak compactness and moment convergence theorems (see e.g. [H]), 
it follows that, for k e J 7 " and all < s < 4 







jk,nPn 



Thus, the term coming in Lyapounov's condition (13.31) is equal to 

n 



/3 (4) 



V 



m n + 



o 
o 



a-p°j 2 

,o ^4, 



/ 



(l-p°) 4 m ra + p°e ra 
[(l-p°) 2 m n + p°£ n ]' 
1 

(l-p°)m n + p°f„ 
1 



1 -p°)m„, 

which tends to as n — > oo. The other error terms given in Theorem 13.21 are 
estimated in the same way. 

The only thing left to check is that oft = 0{ln)' ^ n ^ ac ^' smce In = 
fl(m n ), this relation will only hold when p^ is uniformly bounded away from 1. 
However, for any k G J-° and for any 9 e [— 7r, 7r], 



IVfc,n( )l 



< 



< 




1 + 0(9) 



0(9) 



provided that a < p° n < 1, where a is some fixed constant. This bound can 
be used to replace Equation (IA.2j) in the proof of Proposition 13.11 by 



/ 

J 5, 



-Wx 



: <p n (9)de 



<5„<|6»|<7T 





1 


<- / 


I 1 T 6 . 


J\e\>8 n 



d0 



which is exponentially small in £ n , since 5 n a n = f2(l). 
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5 Towards more tangible assumptions 



The assumptions used in the results of the previous section may seem difficult 
to check in practice. However, as shown hereafter, they can be replaced (at 
the expense of a loss in generality) by simpler properties directly related to 
the service mechanisms of the queues. 

The next lemma provides a realistic context in which IA3I is satisfied. 

Lemma 5.1 Assume that 
(i) there exist sequences R(q) and T(q) such that 

lim ^/R{l)---R{q) = 1, 

q— >oo 

lim T(q) = oo, 

q— >oo 

and, for any q > 0, 

fj-k, n (q) > R{q)vk,m for k e T n , 

y-k,n{q) > T(q)fj, kin , for k E Z n ; 



[ii) there exists a constant B < oo such that 

l0 



A, 



< B, for all k e X n . 



rfeen lA3l holds. 



Remark This lemma can be applied in particular to any mixing of M/M /oo 
and multiple-server queues with at most c servers, with 



R(q) = min 



T{q) = q. 



Proof For each queue k e T n such that p k , n — A and for all r e N, we 
have 

M <-£ q ' A " 



< OO. 



In particular, f k , n (K^k,n) = and 



A„7T fe , 



(A n 7Tfc )r 



An^fc.n 

n^k,n 
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Similarly, for any r G N, 



The same computations can be applied to k G X n , thus proving [A3]- (ii) . 

■ 

The results of Section H] can be easily generalized to a situation where 
some M/M/oo queues of I n become saturated, in which case IA3+ ( ii) is no 
longer satisfied. Indeed, the characteristic function of the number of clients 
X in an M/M/oo queue with parameter p can be written as 



which means that a saturated infinite server queue can be replaced by several 
non-saturated infinite-server queues without changing the distribution of S n . 
Therefore, the results of Section H] still hold, except for marginal distributions 
containing one of the saturated queues. 

Theorems 14.31 and 14.41 also required assumption IA1I on the service mech- 
anisms of the so-called "saturable" queues. It is often enough to restrict 
ourselves to the following two categories of queues, which encompass the 
standard M/M/c queue. 

Lemma 5.2 Assume that, for any k e T^, either 
(i) there is a constant q c , independent of k and n, such that 





(5.1) 



or 



(ii) for some finite constants £ 



mm 



and £ 



,max; 



log 



+ A k , n (q) 



(5.2) 



t^k,n 



with 





max > 



uniformly in k and n. (See also Section^. 
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r/ien[Al] holds. 

Proof In view of Equation 
estimate is related to 



for any fixed k and n, the quantity to 



f (\ nr J>\ V* (An7T fc , n e ie ) g 

/fc,n(.A n 7Tfe in e J = > 7-T 

~ Wfc.n 1 • • • 



g=0 
oo q 



g=0 p= 



J A*fc,n(P) 



For the sake of brevity, let us omit the k and n subscripts and define, for 
any z£C, |z| < 1, 



oo q 



Thus, we have to estimate g(pe l9 )/ g(p), for 6> G [— 7r,7r] and p < 1. This 
proof proceeds in steps: 

a) Assume first that (15.11) holds. Then 
1 



9K Z ) 



l-z 



<jc q q c 

I 1 „q , ^qc + 1 TT ^ 



°( 1 -*)En-7V'+* to+1 IT 



g=U p= 



and Assumption IA1I holds with £ = 1. 
b) Under (I5.2p . one obtains, for q > 1, 



TT — 

y /x(p) 



p =i " p =i 

exp[(£-l)C-A].g^ 1 [l + ^] 

clef \ oo 



where C is the Euler constant, A = J^Li A(p), anc ^ a g * s uniformly 
bounded. In the remainder of the proof, let 

K d = f exp[(e-l)C-A]. 



c) Let, for \z\ < 1 and s G C, 



<f)(z, s) = 
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Then, for Re(s) > 0, 



[Z,S) 



T(s) 



00 t s ~ l dt 



& — z 



In fact, this integral representation can be used to get an analytic 
continuation with respect to s, by introducing the (classical) Hankel's 
contour. This yields, for all \z\ < 1 and Re(s) > 0, 



[z,s) 



ir(i - s ) r (-ty-Ht 



2tt 



e L — z 



Distorting £ to include the zeros of e'— z, the following expression holds, 
for Re(s) < and all values of z such that | arg(— log 2; + 2imi)\ < 71: 

<f>(z,s) = T(l - s)^(-log2 + 2m 7 r) s - 1 . 



d) Using this expression, simple computations yield, when £ > 1 and 
\z\ < 1 



9{z) 



and, finally, 



l + K 

K 
log^ z 



3=1 



K 



+ > + £ 



log z 



log z — 2in7r 



+ log^ 2; 

9=1 



ff(pe ie ) 

q{p) 



1 - pe ie 



it 



[l+£p(e ie -l)]. 



This concludes the proof of the lemma. 



6 Applications 

6.1 A Jackson network with convergence properties 

Consider the basic Jackson network (consisting of M/M/l queues with con- 
stant service rates) analyzed in [10J. 
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In this case, 

m n (tXJ = —— , with r fcjn = . 

Under the assumption made in [10J that the counting measure 



I n (A) = -Card(fc : r fe>n G A), 
n 



defined for all Borel sets A, converges weakly to a probability measure /, we 
have 

m n (tX° n ) f 1 tr 



and 



lim "-"v-n/ = / dI , r \ 

n->oo n In 1 — tr 



lim / dl(r) = X cr < oo. 

t-+i- J 1 - tr 



Thus, the results of [10] are contained in the theorems of Section HI taking 
m° = n\ cr , which is then a strongly critical sequence for C n . 



6.2 A network with tight bottlenecks 

As pointed out in the introduction, there are cases of interest with m n = o(n). 
This will be illustrated in the next example. 

Consider a closed network consisting of s n subnetworks of M/M/l queues 
having each a unique entry point, in which a fixed number m of tasks circu- 
late. The queues are subject to failures, taking place with some probability 
/ < 1. When a failure occurs, the task returns to the entry point of its 
current subnetwork. Tasks visit the various subnetworks according to some 
probability matrix. 

This model exhibits tight bottlenecks, when the number and the size of 
the subnetworks grow. This fact, for the sake of simplicity, will be illustrated 
on a very simple topology, presented in Figure [TJ all subnetworks are associ- 
ated in tandem, and each of them consists itself of l n queues in tandem, with 
unit processing rates. 

Here, the invariant measure of the routing matrix has the form 

where 71"/^ is the invariant probability associated to the fc-th queue of an ar- 
bitrary subnetwork. A straightforward computation, using symmetry prop- 
erties, yields, for any t e]0, 1[, 
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subnet 1 subnet 2 subnet s. 



f 


/ 
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f 


i i 
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/ 


/ 
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f 


/ 


/ 


/ 
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» H 




1 H 


(-4- 



Figure 1: a compound network of tandem queues 



in, -\jfc-l 



m n (tX° n ) = s n J2 1 _ t ( 1 _fik-i - 



Choosing some fixed ug]0,1[ and assuming that £ n — > oo as n — > oo, we 
have 

Um m n (t\° n ) = Lf® 
n^oo m n {u\° n ) Lfiu)'' 

where Lf is defined on ]0, 1[ as 



and lim^i_ = oo. 

Therefore, m n (-uA^) is a strongly critical sequence for the network and 
the size of the queues remain uniformly bounded if, and only if, 

m n = (m n (uX° n )) = 0(s n ) = o(n). 



6.3 A service vehicle network 

Consider a fleet of vehicles serving an area consisting of n stations forming 
a fully connected graph. These vehicles are used to transport goods or pas- 
sengers. Vehicles wait at stations until they receive a request, in which case 
they go to an other station. The routing among stations is done according 
to some routing matrix P n . When a request arrives to an empty station, it is 
immediately lost. The request arrivals form a Poisson stream at each queue. 

We model this system as follows: for all < k < n, station k is repre- 
sented as a single-server queue with service rate jik,n which is equal to the 
arrival rate at station k, since arrivals are lost when the station is empty. 
When a vehicle leaves station k, it chooses its destination according to the 
Markovian routing matrix P n = (pke, n )- The duration of the journey be- 
tween two stations k and I is represented by an infinite server queue placed 
on the edge between them. The service rate of this queue when there are 
q vehicles traveling between k and £ is qfike,n- Note that, contrary to the 
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convention used throughout this paper, the total number of queues is n 2 + n. 
Let (vr l n , . . . , 7r n ri ) be the invariant measure of P n , defined as in (12.1 1) . Then, 
with obvious notation, for all k, I G [1, n], for all 9 G [— 7r, 7r], 

dcf ^n^k,n dcf ^n^k,nPkl,n 

Pk,n — "7j > PH,i 



dcf Pfc,n dcf 

^fc,n — ^ , m H,n — Pfc£,n, 

1 — Pk,n 

(~\ n \o~' lrn k n@ 

1 - PM e 

Define as in SectionHJand assume that its cardinal is some fixed integer 
K > 1. Lemmas 15.11 and 15.21 apply, taking -R(g) = 1, T(q) = q and £fc jn = 1 
for q > 1 and k G J 7 ". Thus, when lA4l holds. Theorem 14.31 can be used and 
estimates of many performance measures can be derived, with corresponding 
error terms. 

Some questions of interest arise: 

• which maximal efficiency can be expected from this system? 

• how many vehicles should be provided? 

To answer these questions, it is convenient to define the loss probability 
as n 



^loss(^) 



Era 
k=l Pk,n 

'Pioss(^) is the proportion of customers that are lost because they arrive at 
an empty station. This is a good indicator of the quality of service provided 
by the network. Under appropriate conditions as n — > oo: 

121=1 Vh,nH X k,n = °) 



Vi oss (n) 



Era 
k=l Pk,n 

A " (6-1) 



2 X)fc=l Pk,n 



The last expression is a decreasing function of A n , which is itself bounded 
by A° . Therefore, the minimum loss probability is attained when A n — > A° ; 
this happens with 



m n = (1 - 9 n )rh° n , lim 9 n = 0, 

n— >oo 



where 9 n is chosen to satisfy the assumptions of Theorem I4.31 (z) . With this 
choice of m n , (16. ip holds with 

X n = X° n (l + 0(9 n )), 
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which is asymptotically optimal. Consequently, a "good" value for m n is 
m n = rh° n , and having a number of vehicle proportional to the number of 
stations can be a poor choice, especially when some stations are more loaded 
than others. These stations act as bottlenecks of the system, which should 
be removed by altering the routing probabilities. 

7 General remarks 

First, a chief difficulty of the analysis is due to the need of dealing with 
rate of convergence and limits of densities: this is the field of Berry-Esseen 
theorems and large deviations. 

Secondly, the results have been obtained under several technical assump- 
tions (especially uniformity), which in some sense are unavoidable. This 
means precisely that the choice of conditions slightly different from IA1[ IA3I 
and IA4I would have led to different families of limit laws having infinitely 
divisible distributions. 

In particular, from a physical point of view, it is worth commenting on 
equation (15. 2p . The inequality ^ )Tl > 1 implies that the maximum service rate 
of the queues in T° n is reached from below; this is not the case if < < 1? 
and the analysis was omitted, since the technicalities involved would have 
made the text unnecessarily obscure. At last, the case < dealing with 
other types of singularities (for instance logarithmic), was not carried out, 
and would yield other limit laws. 

The future class of problems of interest concerns some non-product form 
networks. 

A Appendix 

A.l A bound on periodic characteristic functions 

One of the problems arising in the computation of convergence rates in the 
Central Limit Theorem is to find upper bounds on the modulus of a charac- 
teristic function <p(6) for 9 away from 0. One typical property used can be 
stated as follows: 

there exist 9q > and a < 1 such that, for all \9\ > 9q, \<p{9)\ < a - 

It is pointed out in Feller [5] that this condition is usually easy to fulfill 
in practice, as long as X does not have a lattice distribution. Unfortunately, 
we are in the lattice case and thus must cope with the periodicity of <p. 

Next lemma shows how a bound on \<p(9)\ can be derived for \9\ < n. 
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Lemma A.l Let X be an integer-valued random variable with distribution 
P{X = k) = p k , k E N. Define 



2 dcf \ ^ 

7 =2^ 



k=0 



P2kP2k+l , . / 1\ 

< mm Var A, - I , 

P2k+P2k+1 V 4/ 



where the summands are taken to be zero when p2k = P2k+i = 0. Then, for 
any 9 6 [— 7r, 7r], i/ie characteristic function tp of X satisfies: 



Proof We have 



\<P 



Moreover, 



P2k + P2k+ie 



Hence, for 6 E [0, 7r] 



|^)|<exp(-^ 2 ). 



fc=0 



fc=0 



(P2k + P2k+i cos 6 1 ) 2 + p 2 2k+l sin 2 

a/ {P2k + P2k+l) 2 - 2p 2 fcP2fc+l(l - COS( 
P2kP2k+l 



< P2k + P2k+1 ~ 



P2k + P2k+1 



(1 — COSI 



(A.1) 



< i-(i-cos0)E 



fc=0 



P2kP2k+l 
P2k + P2k+l 



2 
~ 2 



< exp( 



7rf_ 

7T 2 



e 2 



which yields (lA.lj) . That 7 2 < VarX can be seen by a Taylor expansion of 
if in the neighborhood of 9 = 0, while the relation r j 2 < 1/4 follows from the 
trivial inequality 



P2kP2k+l P2k + P2k+l 



P2k + P2k+1 



7 has the desirable property to be zero when X is an integer variable with 
a span strictly greater than 1, in which case the period of if is less than 2tt. 
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Another desirable property would be that 7—7-00 when the moments of X 
are unbounded; since 7 < 1/2, this is obviously not possible here. That this 
"feature" is somehow unavoidable can be seen on the following example: 

2 + e id lA e ike 



1 - e ld ) ln(l - e w ) 



1 + e 



2 

The random variable having ip as characteristic function admits no finite 
moment of order greater or equal to 1, but no bound on \ip\ is substantially 
better than fjA.lj) . 



A. 2 Proof of Propositions 13.11 and I3T3 



Proof of Proposition 13.11 Using a Fourier inversion formula, the left hand 
side of (13. ip can be rewritten as 



A6x 



2tt J_ w r ' iy ' 2tt 
Thus, our goal is to evaluate the quantity 

r-5,. 



e '"""e 2 du. 



00 cp-0 1 



-Wx 



<Pn{0) ~ e 



00 

T 2 S 2 

2 



CT 2 e 2 

e - iex e -^d9 + 



It is known that 



|0|><s„ 



e 2 d9 « 



t20 



|6[<5„,ir] 



2 .2 



e 2 



applying Lemma LA. II to ip n , we get 



-Wx 



Vn{9)d9 



s n <\e\<TT 



< 



\8\>S n 



< -^dH ^ ( ;( J-e" 2 ^ 



Finally, we obtain a bound on \l n \ which is uniform in x: 

«8 n 



< 



e 2 



+ 



2 <;2 



e 2 



dO 



+ 



-e 5 



(A.2) 



(A.3) 
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We proceed now to estimate the above integral, so that implicitly \9\ < 5 n . 
The derivation relies on the following simple inequality, valid for all complex 
numbers X\, . . . ,x n and yi, . . . , y n \ 



Fi • • -x n - yx ■ ■ -y Tl 



< 



k=l 



xi ■ ■ -Xk-iWxk - yk\\yk+i ■ ■ 'Vr, 



(A.4) 



which will be used with Xk = fk,n(9) and yu = exp(—al n 9 2 /2). 

The characteristic function tpk,n of the random variable Xfc >n satisfies (see 
for example Loeve [H]) 



n2 



(2+r) 



— Pk,n 



\e\ 2+r 



(A.5) 



Hence, using the inequality \e x — 1 + x\ < x s /s, valid for all x > and 
1 < s < 2, 



— lH 



< 



2 



< A 



<Pk,n(0) ~ 1 + °k,n — 
01 2+r 



+ 



(2 + r) 

k.n 



a 



2+r 
k,n 



2 

e\ 2+ ' 



& 2 



2 - 1 + <n~ 



< C r) l^l 2+r - 



2 

(A.6) 



To find an upper bound for |<£>fc, n |, assume first ak,n&n < 1, so that 

\e\ 2+r 



n2 



(A.7) 



In fact, (IA.7P also holds when a^, n b n > 1, since in this case 

-<n + * ~<n + VIZ*: > 0- 

From (13 .4p . we can choose n such that a^n < a n /2 and, using (IA.4p . (IA.6[) 



and (IA.7I) . we find 



7 2 e 2 

J n u 
1 



fe=l 

< ^ r) \0\ 2+r ^p(-^j)- (A.8) 



Equation (13. ip follows, since the integral in (1A.3|) is bounded by 



Sn 



<p n (9) - e 2 



d9 < (3 n 2+r) j |#| 2+r exp(-CT, 

J — oo 



o 



1 



oo 

(2+r) 



d0 



2+r 
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The proof of (13.21) of the proposition is similar, although the computations 
be more involved. Redefine /„ as 



I„ = I e ~ Wx Vn (6)d6 



e- Wx (l-iP ( n 3) -)e-^de, 



To find a bound for \I n \, we have to estimate 

- 9' c 



(A.9) 



< 



<p n (6) - e-—~ l 



2 /)2 a (3) fl 3 



1 + 



e 2 



The first part of the r.h.s. of f]A.9|) is evaluated as above with (IA.4|) 
and f]A.7|) replaced by 



^,nW<exp(-a fc 2 n + /3^ ri ) 



2 



For the second part, we use the following inequality, valid for r > (see 
e.g. Loeve [H]) 



r «(3) 1 1+ 



< 



B (:: 



(3+r) 



a 



3+r 



which yields 



n3 

1 + 



< 



fn 6 



1+5 /3(3+'') n 3 + r \0\ 3 + r 



< 



6 



and ([32]) follows. ■ 

Proof of Proposition 13.31 The proof of this proposition is similar to the 
proof of Proposition 13.11 and is only sketched here. Define 

/ x def 1 
LO(U) 



1 — iu ' 

do_f Sje-FO 171 j,n + x 

J 

CXn. 



and 



In = a n r e- Wa ^cot(9) J] ip k>n (8)<p n (6)d9 
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u 



+ 
+ 



-i«2M, ,£ 



-7ra r , 



ul n (u/a n ) - w €n (u) 



du 



e °% 2 du 



(A.10) 



|li|>7TQ n 



The evaluation of these integrals depends on the following straightforward 
estimations, valid for \u\ < 7ra„, 



\ut(u/a n )\ = O 
\u^{u/a n )-uj^{u)\ = O 



O 



(1 + U 2 )^l 2 

1 n 2 
a n (1 + u 2 )€™ 
1 + |u 



and on flA.8[) . which yields for \u\ < a n 5 n , 



(p n {u/a n ) - e °>n 



O 



a 



(p n {u/a n ) 



I X] \ a 2 8 



ol u 2 



Moreover, we use the following approximation, valid for a, b > and for 
sufficiently small z: 



J(a,b,z] 



Z" 00 17/1 
dcf / I U I 



Loo (l + « 2 

These relations, together with OA.lOp . yield: 



e~ z u du = 0(l) + 0(z 2b - a - 1 ). 



+ 



+ 



—)j(2,Z n 



CT r . 



Or. 
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+ 



a,, 



v sn e 2 av 



V>&n% 



Oil- + ^H(^) 2+r + Pn +r) (*n\^ 



a,, 



+ ( 



a 2 + r \a z 
e 



2+r 



a, 



f n <Jn U-n 



To conclude the proof of (13. 8p . the second term coming in the definition 
of l n is evaluated using Parseval's identity and classical tools of complex 
analysis (see e.g. Lavrentiev and Chabat [8]). This yields 



e~ lJ/nU u* n (u)e °£ 2 du 



y&i lg Vn 



l + 0(% 
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